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ABSTRACT 

The fine structure of solar prominences and filaments appears as thin and long 
threads in high-resolution images. In Hor observations of filaments, some threads can 
be observed for only 5-20 minutes before they seem to fade and eventually disap- 
pear, suggesting that these threads may have very short lifetimes. The presence of an 
instability might be the cause of this quick disappearance. Here, we study the thermal 
instability of prominence threads as an explanation of their sudden disappearance from 
Ha observations. We model a prominence thread as a magnetic tube with prominence 
conditions embedded in a coronal environment. We assume a variation of the physical 
properties in the transverse direction, so that the temperature and density continuously 
change from internal to external values in an inhomogeneous transitional layer repre- 
senting the particular prominence-corona transition region (PCTR) of the thread. We 
use the nonadiabatic and resistive magnetohydrodynamic equations, which include 
terms due to thermal conduction parallel and perpendicular to the magnetic field, ra- 
diative losses, heating, and magnetic diff"usion. We combine both analytical and nu- 
merical methods to study linear perturbations from the equilibrium state, focusing on 
unstable thermal solutions. We find that thermal modes are unstable in the PCTR for 
temperatures higher than 80,000 K, approximately. These modes are related to tem- 
perature disturbances that can lead to changes in the equilibrium due to rapid plasma 
heating or cooling. For typical prominence parameters, the instability time scale is of 
the order of a few minutes and is independent of the form of the temperature profile 
within the PCTR of the thread. This result indicates that thermal instability may play 
an important role for the short lifetimes of threads in the observations. 



Subject headings: Sun: filaments, prominences — Sun: corona — Instabilities — 
Magnetohydrodynamics (MHD) — Magnetic fields 
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1. INTRODUCTION 



Solar prominences and filaments are large-scale magnetic structures of the solar corona. The 
main issues regarding the physics, dynamics , and modeling of t h ese c oronal inhabitants have been 
recently reviewed by Labrosse et al.l (120101) and iMackay et al.l (|201Q) . High-resolution observa- 
tions reveal that prominences and filaments are formed by long (5" - 20") and thin (0".2 - 0".6) 
fine structures, usually called thread s. Although the existence of the fine struct ure of prominences 
was discovered long time ago (e.g., iMenzel & Wolbachlll960l: lEngvoldlll976|) . its properties and 
dynamics could only be studied in more detail with recent high-resolution observations. The fine 
structures show up as dark ribbons in Ha images of fil a ments on the solar disk from the Swedish 
Solar Telescope (e.g., lLinll2004l: Lin et allboovi boosi bo09L and as bright features in observa- 
tions of prominences in the solar limb from the Solar Optical Te l escope aboard the Hinode satellite 



(e.g.. lOkamoto et al.ll2007l : iBerger et al.ll2008l: IChae et al.ll2008l : iNing et al.ll2009l: ISchmieder et al 



2010h . Statistical studies show that the orientation o f threads w ith respect to the filament long 
axis can significantly vary within the same filament (ILinl 12004 1). with 20 degrees a mean value 
typically reported. Vertical threads are more commonly seen in quiescent prominences (e.g.. 



Berger et al.l 120081: IChae et al 



2008h whereas h orizontal threads are ii sually observed in active 



region prominences (e.g., lOkamoto et al.ll2007|) . ISchmieder et al.l (|2010() recently pointed out that 
vertical threads might actually be a pile up of horizontal threads which seem to be vertical struc- 
tures when projected on the plane of the sky. However, this is still a matter of controversy. Since 
threads are observed in both s pines and barbs, it is believed that they are the basic building blocks 
of prominences and filaments (^Engvoldll2004^) . 



Theoretically, the fine structures have been modeled as magn etic flux tubes anchored in the 
solar photosphere (e.g., Ballester & Priest|[l989 : Rempel et al. Ill999l) . which are only partially filled 
with the cool (~ 10^ K) filament material, while the rest of the tube is occupied by hot coronal 
plasma. Therefore, the magnetic field is oriented along the axis of the fine structure. This model 
is conceptually in agreement with the idea that the dense prominence material is trapped in dips 
near the apex of a magnetic arcade connecting two photospheric regions of opposite magnetic 
polarity. The dips are supposed to correspond to the observed threads, which are piled up to 
form the prominence body. It has also been suggested from differential emission measure studies 
that each thread might be surrounded by its own prominence-corona transition region (PCTR) 
where the plasma phys ical properties would abruptly vary from prominence to coronal conditions 
dcirigliano et al.lboo4 . 



Prominence threads are highly dynamic (see, e.g.. lHeinzelll2007l : lEngvoldll2008|) . For exam- 
ple, transverse thread oscillations and propagating waves along the threads seem to be ubiquitous 
in prominences, wh ich have been interpreted in t erms of magnetohydrodyn amic (MHD) waves 
(see the reviews by iBallesterl 120061: lOliverl 120091: lArregui & Ballester! 1201 0|) . Mass flows along 
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threa ds, with typical flow ve l ocities of less t han 30 km s \ h ave been also frequently reported 
(e.g. Jzirker et alJll994[l998HLin et"aDl2003Hchae et aljboosl) . An interesting property of the ob- 
serva tions is the apparent short l i fetime of some threads when they are observed in Hor sequences 
(e.g., Einll2004l : Ein et all llOOsi l2009l) . Typically, the threads can be followed for only 5-20 
minutes before they seem to fade with time and eventually disappear. The cause of this quick dis- 
appearance is unknown, although several explanations hav e been proposed. A possible explanation 



is rela ted to the presence of flows and mass motions (e.g. JChae et al.ll2008U2010l: 



Schmieder et al. 



20101) . which may trigger a Kelvin-Helmholtz instability (KHI). Recently. ISoler et al.l (l2010bh in- 
vestigated t he KHI in magnetic flux tu bes due to shear flows generated by transverse motions of the 
tube, while IZaqarashvili et al.l (|2010l) studied the KHI in twisted tubes due to longitudinal flows. 
When the results of both papers are applied to prominences, one obtains that neither the observed 
transverse velocity amplitudes of oscillating threads nor the longitudinal flow velocities are large 
enough to trigger a KHI on time scales consistent with the apparent l astin g time of the thr eads. 
On the other hand, a different explanation was suggested by Lin ( 2004 ^ and Lin et al. ( 2005 ). Ac- 
cording to these authors, rapid cooling or heating may cause the maximum of the plasma emission 
to fall outside the bandpass of the filter, and so the thread would become invisible in Ha images. 
A mechanism that could lead to a rapid heating or cooling of the prominence material is a ther- 
mal instability. Here, we explore this possibility by studying the thermal instability of prominence 
threads 

Thermal or condensation modes have bee n extensively investigated in homogeneous plasmas 
(e.g., Parker|[l953 : Fieldlll965 ; Hevvaertsll974 ). A relevant work in the context of prominences was 
performed by lCarbonell et al.l(|2004[) . who studied the thermal mode in homogeneous plasmas with 
prominence, PCTR, and coronal conditions, considering parallel thermal conduction to rn a gnetic 
field lines and the optically thin radiative loss function of iHildneii (|1974l) . ICarbonell et al.l (|2004|) 
obtained that, for long wavelengths, the thermal mode is unstable for PCTR temperatures since 
thermal conduction is not efficient enough to stabilize the thermal disturbanc e. In the case of an in 



homogeneous plasma, thermal mode s were studi ed in detail in th e work s by IVan der Linden et al 



(1199 ih . IVan der Linden & GoossensI (|l99l|) and IVan der LindenI (| 19931) . As the present investi- 
gation is based and inspired by these previous works, we summarize their relevant results in the 
following paragraph. This will also enable us to clearly define the advancement made in the present 
paper compared with these previous investigations. 



Van der Linden et al.l (|l99lh studied the linear nonadiabatic MHD spectrum of a magnetic 



cylinder taking into account the effect of plasma inhomogeneity in the transverse direction. They 
included radiative losses, heating, and thermal conduction parallel to the magnetic field lines but 
neglected perpendicular thermal conduction. These authors pointed out the important result that, 
along with the classical Alfven and slow continua, there exists an additional thermal continuum, 
which can be unstable depending on the physical properties of the equilibrium. Subsequently, 
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Van der Linden & GoossensI (|l99l|) investigated the effect on the thermal continuum of themial 
conduction perpendicular to the magnetic field lines, and showed that perpendicular thermal con- 
duction replaces the continuum by a dense set of discrete quasi-continuum modes. These quasi- 
continuum modes retain the basic stability properties and growth rates of the continuum, and are 
confined within thin conductive layers around the position of the continuum singularities. The 
width of the conductive layer, and so the spatial scale of the quasi-continuum mod e s, dep ends 
on the value of the pe rpendicular thermal conductivity. IVan der Linden & GoossensI (|l99lh and 
Van der LindenI (Il993h applied these results to the context of solar prominences and showed that 
the spatial scales related to the most unstable quasi-continuum modes of the spectrum are con- 
sistent with the size of the prominence threads reported from high-resolution observations. This 
suggests that perpendicular thermal cond uction may be responsible for the formation of promi- 
nence fine structure. In subsequent works. Ilreland et al.l(|l992lll998|) investigated the influence of 
magnetic diffusion on the thermal continuum. These authors concluded that magnetic diffusion 
plays a similar role to that of perpendicular thermal conduction as both mechanisms replace the 
continuum by discrete modes. 



In the present work we follow the method of I Van der Linden et al.l (|199lh and lVan der Linden & Goossens 



(|1991i) and investigate the thermal stability of an inhomogeneous prominence thread model. The 
equilibrium model adopted here is a cylindrical magnetic tube, representing a prominence thread, 
embedded in a coronal environment. The plasma physical properties, i.e., temperature, density, 
etc., are inhomogeneous in the transverse direction, and vary from prominence to PCTR and coro- 
nal values. In this investigation, we assume that the equilibrium is uniform in the longitudinal 
direction. The effect of longitudinal plasma inhomogeneity is relegated to a forthcoming study. 
We use the nonadiabatic and resistive MHD equations and superimpose linear perturbations on 
the equilibrium state. Parallel and perpendicular thermal conduction, radiative losses, heating, and 
magnetic diffusion are the nonideal effects included in our equations. We combine both analytical 
methods and numerical computations to study the properties of the discrete thermal continuum 
modes, focusing on the unstable modes with largest growth rates. 

This paper is organized as follows. Section |2] contains the basic equations and a description 
of the equilibrium configuration. Then, Sections |3]-|5]contain a theoretical study of the properties 
of the thermal instability. The stability of the thermal continuum in the absence of perpendicular 
thermal conduction and ma gnetic diffusion is discussed in Section [2]by following the analysis of 
Van der Linden et al.l (1199 ih . while the unstable thermal modes are investigated in Section |4] when 
both perpendicular thermal conduction and magnetic diffusion are included. A parametric study 
of the growth rate of the most unstable solution is performed in Section \5\ Then, we discuss 
in Section [6] the physical implication of our theoretical results for the stability and lifetime of 
prominence threads. Finally, the summary of our results is given in Section|7l 
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2. BASIC EQUATIONS AND EQUILIBRIUM 

The basic MHD equations governing the dynamics of a nonadiabatic and resistive plasma are 

§^.pV.v = 0, (1) 

p^ = -Vp + -(VxB)xB, (2) 
Dt fi 

^ = V X (v X B) - V X (77V X B) , (3) 
ot 

^ - y ^ + (r - 1) [pL(T,p) -V-iK-VT)- 77IV X Bf] = 0, (4) 

P = pR^, (5) 

along with the condition V ■ B = 0, where ^ = ^ + v • V is the material derivative for time 
variations following the motion. In Equations (d])-©, p, p, T, v, and B are the mass density, gas 
pressure, temperature, velocity vector, and magnetic field vector, respectively. In addition, p. = 4nx 
10"^ N A"^ is the magnetic permeability, y = 5/3 is the adiabatic index, R = S.3 x 10^ m^ s"^ K"^ 
is the ideal gas constant, L{T,p) is the heat-loss function, k is the thermal conductivity tensor, 77 is 
the magnetic diffusivity, and jx is the mean atomic weight. In a fully ionized medium, fi = 0.5. 



2.1. Nonideal terms 

The induction (Equation (O) and energy (Equation dD) equations contain several nonideal 
terms whose physical meaning is explained next. The term with the factor 77 in Equation ([3]) 
corresponds to Ohm's magnetic diffusion, while the equivalent term in Equation (jH) accounts for 
Ohm's heating. The other nonideal terms in Equation (H) correspond to nonadiabatic effects, i.e., 
thermal conduction, radiative losses, and an arbitrary heating input. 

As for thermal conduction, the thermal conductivity tensor, k, can be expressed in terms 
of its parallel, k\\, and perpendicular, k±, components to the magnetic field, namely k = KwegeB + 
/Cx (l - eBes), where eg = B/ |B| is the unit vector in the magnetic field direction and I is the identity 
tensor. The parallel and perpendicular conductivities are dominated by the effect of electrons and 
ions, respectively. The perpendicular electron conductivity and the parallel ion conductivity are 
negligible in fully ionized plasmas. Expressions for the conductivities in MKS units for a hydrogen 



-6- 



plasma are (see, e.g.. |Parkenll953l: ISpitzerlll962t [Braginskiilll965|) . 

T5/2 



'<\\ 



1.8 X 10 



-10' 



= 1.48x10 



^42 



In A 
lnAp2 



(6) 
(7) 



where nii is the proton mass and In A is the Coulomb logarithm, whose val ue is general ly between 
5 and 20 and has a weak dependence on temperature and density (see, e.g.. |Priestill982|) . 



Finally, the heat- loss function L(T,p) accounts for the balance between radiative cooling and 
heating. The determination of an analytical function of the temperature and density that describes 
radiative losses of the prominence plasma is a very difficult work that requires the numerical solu- 
tion of nonlocal thermodynamic equilibrium (NLTE) radiative transfer equations. This is beyond 
the purpose and scope of the present study. One reasonab le semi-empiric al approximation to an 
expression for the radiative los s function was obtai ned by iHildnerl (|l974() . This author assumed 
an optically thin plasma (e.g., ICox & Tucker] 1 1969|) and performed a piecewise fit of the radia- 
tive losses previously computed by several authors as a function of ter nperature. Here , we adopt 
Hildner's approach. The functional expression of L{T,p) considered by lHildnei](|l974|) is 



L{T,p)=px*T"-h, 



(8) 



where x* and a are piecewise constants depending on the temperature, and h is an arbitrary heat- 
ing function. The assumption of an optically thin plasma seems a reasonable approximation for 
coronal and prominence-corona transition region (PCTR) tempera tures, whereas cool prominence 



plasm as may be considered optically thick. Some authors (e.g.. iRosner et al.lll978 



Milne et al. 



119790 have proposed corrections to the values of x* and a in the range of cool prominence tem- 
peratures, i.e., T < 15,000 K, in order to represent radiation losses in optically thick plasmas 
using Equation ([8]). In the literature, there are more recent parametrizations for x* and a which 
update Hildner's values. In particular, the so-called Klimchuk-Raymond parametrization (see, e.g., 
Klimchuk & CargilllllOOlh may be a better representation of the radiative losses in the hotter part 
of the PCTR and in the solar corona. On the contrary, the Klimchuk-Raymond function may not be 
adequate in the cool part of the thread. We use the Klimchuk-Raymond function as an alternative 
to Hildner's parametrization in the hotter part of the equilibrium. The values of the parameters x* 
and a for various temperature ranges and regimes are given in Tabled] 



2.2. Equilibrium configuration 



Equations ([I])-© are applied to the following equilibrium configuration. The prominence 
thread model is composed of a straight and cylindrical magnetic flux tube embedded in a coronal 
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Regime 


Temperature range 




X* 


a 


Prominence- 1 . 1 


T < 15x10^ K 




1.76 X 10"^^ 


lA 


Prominence- 1.2 


T < 15x10^ K 




1.76 X 10"" 


11 A 


Prominence- 1.3 


T < 15x10^ K 




7.01 X 10-^°"* 


30 


PCTR-2 


15x10"^ K < r < 8x10^ 


K 


4.29 X 10^" 


1.8 


PCTR-3 


8x10^ K< T <3xlO^ 


K 


2.86 X 10^'' 


0.0 


PCTR-4 


3x10^ K< T < 8x10^ 


K 


1.41 X 10" 


-2.5 


Corona- 5 


T > 8x10^ K 




1.97 X 1024 


-1.0 


Klimchuk-Raymond- 1 


T < 10^ ''^ K 




3.91 X lO'' 


2.0 


Klimchuk-Ray mond- 2 


J 04.97 ^ 7^ < ]^q5.67 ^ 




3.18 X 10^4 


-1.0 


Klimchuk-Raymond- 3 


T > 10^-^^ K 




6.81 X 10'^ 


0.0 



Table 1 : Values in MKS units of the parameters in the radiative loss function (Equation ([8])) cor- 
responding to several temperature regimes. Pr ominence- 1.1, PCTR- 1 , PCTR-2, and PCTR-3, and 
Corona-5 regimes are para metrizations from Hildnerl (119741) . Promine nce- 1 .2 and Prominence- 
1.3 regimes are taken from IMilne et al.l (|l979h and iRosner et al.l (|1978b . respectively. The three 
Prominence regimes represent different plasma optical thicknesses. Prominence- 1 . 1 correspond- 
ing to optically thin plasma, while Prominence- 1.2 and Prominence- 1.3 are for optically thick 
and very thick plasmas, re spectively. The three Klimchuk-Raymond regimes are adapted from 



Klimchuk & Cargilll (120011) . where we only have taken into account the range of temperatures 



considered in our equilibrium. 



environment. We use cylindrical coordinates, namely r, ip, and z for the radial, azimuthal, and lon- 
gitudinal coordinates, respectively. In this configuration, the equilibrium quantities are invariant 
in the azimuthal and longitudinal directions, so they depend on the radial direction only. Here- 
after, the equilibrium quantities are denoted by a subscript 0. The equilibrium magnetic field is 
straight and homogeneous, Bq = B^e^, where, for simplicity, Bq is the same constant everywhere. 
Magnetic twist and longitudinal plasma inhomogeneity are effects that will be included in future 
investigations. 

Since the magnetic field is straight and homogeneous, the equilibrium gas pressure, po, is also 
homogeneous. The equilibrium density, po, and temperature, Tq, profiles must verify the condition 
of energy balance, which from Equation (H) is 

where we have assumed that there are no flows in the equilibrium, i.e., Vo = 0. Note that for the 
present equilibrium magnetic field. Ohm's heating term is absent from Equation (|9]). The density is 
related to the temperature and gas pressure through Equation ([5]). We can use Equation ^ to obtain 
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the temperature profile and later use Equation (|5) to compute the density profile. Unfortunately, 
the heating function h is unknown in prominences and in the corona, which makes it impossible to 
apply this procedure. In stead, we follow an alternative method that has been adopted in a number 
of previous works (e.g., IVan der Linden & Goossenslll99ll: llreland et al.ll 199211 1998|) . We choose 
an ad-hoc temperature profile and compute from Equation (|9) the corresponding heating function 
that satisfies energy balance. The temperature prof ile adopted in this wo rk is in agreement with 
theoretical models of prominence threads (see, e.g.. ICirigliano et al.ll2004l Fig. 12). 



Thus, the equilibrium temperature profile is 



To(r) = 



T 

^ P' 

T^PCTR (f ) ■ 

T 



if r<R-ll2, 

if R-l/2<r <R + 1/2, 

if r>R + l/2. 



(10) 



We consider that the prominence thread is composed by a central and cool region with homoge- 
neous temperature, Tp, surrounded by a transverse transitional layer where the temperature abruptly 
increases until the homogeneous c o ronal medium with temperature Tc is reached. Following the 
idea explained by ICirigliano et al.l (|2004|) . this transitional layer represents the particular PCTR 
of the prominence thread. The mean radius of the thread is denoted by R and is located at the 
center of the PCTR. The thickness of the PCTR is denoted by /, and corresponds to the region 
R - 1/2 < r < R + 1/2. So, the cases l/R = and l/R = 2 correspond to a thread without 
transitional layer and a radially full inhomogeneous thread, respectively. In Equation (flOl) and in 
the following expressions, the subscripts p, PCTR, and c denote the central prominence region, 
the transitional layer, and the coronal medium, respectively. Unless otherwise stated, we consider 
the following values of the equilibrium quantities: Tp = 8000 K, = 10^ K, 5o = 10 G, and 
R = 100 km. The central prominence density is fixed to pp = 5 x 10"^^ kg m~^, thus the gas 
pressure is po = 6.64 x 10"^ Pa, which corresponds to /3 = po/ ^5^/2//^ = 0.008. 

In order to assess whether or not the form of the temperature profile influences the results, we 
consider three different expressions of Tpctr (''), which correspond to a linear profile. 



TpcTR (r) = Tp + ^^-j^ (r-R + 1/2) . 



a sinusoidal profile. 



(11) 



(12) 



and a Gaussian profile 



T'pcTR (r) = Tc - (Tc - Tp) exp - 1"— ^y"^ j 



(13) 
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Fig. 1. — (a) Equilibrium temperature profiles, Tq, normalized to the internal temperature, Tp, in 
the case l/R = 1. The different line styles correspond to the linear profile (solid), the sinusoidal 
profile (dashed), and the Gaussian profile (dot-dashed), (b) Corresponding density profiles, po, 
normalized to the internal density, pp. The vertical dotted lines in both panels denote the boundaries 
of the transitional zone. Note that the vertical axis in panel (b) is in logarithmic scale. 

In the case of the Gaussian profile (Equation (fT3l)). note that the exponential factor is exp (-4) x 
0.018 for r = R + 1/2, so the temperature is not strictly T^. when the corona is reached. For 
consistency. Equation (fT3l) also applies for r > R + 1/2 in order to avoid this small jump of the 
temperature when the Gaussian profile is used. Figure H] displays a plot of the three different 
temperature profiles assumed in this work and their corresponding density profiles. 



2.3. Linear perturbations 

We superimpose perturbations on the equilibrium state. Hence, pi, Ti, and pi denote the 
density, temperature, and gas pressure perturbations, respectively, whereas Vi = (vr,v^,v^^ and 

bi = {br,b^,b^ are the velocity and magnetic field perturbations. Next, we assume that these 
perturbations are small, so we restrict ourselves to the linear regime and Equations ([I)-(l5]) are 
linearized. Since the equilibrium is invariant in the azimuthal and longitudinal directions, we 
write all perturbations proportional to exp {st + imi^ - ihz), where m and are the azimuthal and 
longitudinal wavenumbers, respectively, and s is the growth (or damping) rate of the perturbation. 
Equations ©-(EI) become, 

SPl = -poVr - Po \V', + -Vr + — + k,vA , (14) 
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1,1 m \ 

- Bo\v,+ -Vr + -V^j 

(1 iTt^ tn k 
b'' + -K + —kjbip - —br - KK 
r ' r r 

+ i(b[-k,br). 



(20) 



—Tx-{y-\)—pi 
To Po 



- (r-1) 
+ (r-1) 

+ (7-1) 



Po 



po 



Po 



m 



1 



T' 



K^- + <| + K^r^ + (/c,, - k^) ^Kb 



(21) 



where the prime denotes the derivative with respect to r, = — is the sound speed squared 



= is the Alfven speed squared, and the quantities a>p,ojT, and k_i_ are defined as follows. 



PO 



(Op = 



Po 
Po 



L{pq,Tq) + po\ — 



POrj, (dL 



(22) 
(23) 



PoJo 
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Equations (fT4)) - ((2T)) form an eigenvalue problem, with s the eigenvalue. To eliminate one 
variable and simplify matters, we have used the linearized version of Equation (|5]) to write the 
gas pressure perturbation, pi, in terms of the density, pi, and the temperature, Ti, perturbations. 
In addition, to express Equations (fT4)) - (|2T)) in terms of real quantities, we have performed the 
substitutions v^, /Vj. v^, and ib,- — > br. Note that the term related to Ohm's heating in 

Equation (H)) does not contribute in the linear regime and so it is absent from the linearized energy 
equation (Equation (1271)). In addition, since the heating function h represents an external, arbitrary 
heating source, no perturbations of this function have been considered. 

We numerically solve Equations (fT4]) - (|2TI) a nd compute t he eigenvalues and their correspond- 



ing perturbations by means of the PDE2D code (iSewelll 12005 ) based on finite elements. The nu 



merical integration of Equations (fT4])-(l2TI) is performed from the thread axis, r = 0, to the finite 
edge of the numerical domain, r = r^ax, which is located far enough to obtain a good convergence 
of the solution and to avoid numerical errors. We use a nonuniform grid with a large density of 
grid points within the inhomogeneous transitional layer, where the equilibrium properties vary by 
several orders of magnitude. The nonuniform grid also allows us to correctly describe the small 
spatial scales of the eigenfunctions within the transitional layer. The PDE2D code uses a collo- 
cation method and the generalized matrix eigenvalue problem is solved using the shifted inverse 
power method. The output of the program is the closest eigenvalue to an initial provided guess and 
its corresponding perturbations. 

In the present investigation, we restrict ourselves to the study of thermal quasi-continuum 
eigenmodes. For real and m, these solutions correspond to real values of s. When the eigen- 
m ode is a stable, damped solution, s < 0, while unstable therm al modes have 5 > 0. As shown 



by IVan der Linden & GoossensI (|l99l|) and llreland et al.l (|1992|) . some of the basic properties re 



garding the stability of the thermal quasi-continuum modes can be deduced from the properties of 
the thermal continuum in the absence of perpendicular thermal conduction and magnetic diffusion. 
We perform this analysis in the next Section. 



3. THE THERMAL CONTINUUM 



In our equilibrium configuration, the full spectrum of ideal eigenmodes, i.e., when nonideal 
effects are absent, contains soluti ons related to two different continua, namely the Alfven and 
slow (or cusp) continua (see, e.g., ^ Appert et~aL 1974 ). These continua are modified or removed 
if nonideal effects are taken into account. The effect of the different nonideal mechanisms on the 
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continua has been studied by a number of authors, whose relevant results for our investigation are 
briefly summarized next. 

If magnetic diffusion is included but nonadiabatic terms are neglected, both Alfven and slow 
continua are replaced by a set of discrete modes with cornplex s and whose properties depend 



on the value of the difl'usivity (see, e.g., IPoedts & Kerned 119911) . In the presence of radiation 



losses and parallel thermal conduction but for k± = rj = 0, the Alfven continuum remains un- 
affected and the slow continuum is only slightly modified becoming a nonadiabatic slow con- 
tinuum with compl ex s. In addition, a new t hermal continuum is present, which corresponds 
to real values of s (|Van der Linden et al.l Il991|) . When perpendicular thermal conduction is in- 
cluded, the thermal continuum is r eplaced by a dense set of eigenmodes, i.e., a quasi-continuum 
(|Van der Linden & Goossenslll99l|) . Finally, if both perpendicular thermal conduction and mag- 
netic diffusi on are taken into ac count, the three continua are removed and replaced by discrete 
eigenmodes ([Ireland et al.lll992l) . 



To study the properties of the thermal continuum before solving the full problem numerically, 
we fix /fj^ = ?7 = in Equations (fT4l) - (|2TI) . After a lengthy but straightforward process, it is possible 
to combine Equations (fT?l) - (|2TI) to obtain the two coupled, first-order differential equations of 
Van der Linden et al.l (|199ll Equations (2)-(3)), which depend on the coefficient Co defined as 



y 



(25) 



with C, the following third-order polynomial in s. 



Ct 



y-l 



To 



7 I\Pq I Po 
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Van der Linden et al.l (1199 1|) showed that the roots of Co = 0, i.e., the singularities in their Equa- 
tions (2)-(3), correspond to the three continua. Note that Equations (|25l ) and (|26l ) are independent 
of the azimuthal wavenumber m, and so the three continua do not depend on the value of m. For 
a fixed r, the roots of Co = are two purely imaginary solutions given by = -k^vj^, that corre- 
spond to the Alfven continuum, and the solutions of C; = 0. Since C, is a third-order polynomial 
and the nonadiabatic terms are assumed to be small, the roots of C, = are two complex conjugate 
solutions corresponding to the nonadiabatic slow continuum, and a real root corresponding to the 
thermal continuum. It is straightforward to check that in the ideal case k\\ = cOp = ojt = 0, so Ct 

(27) 



becomes 
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Fig. 2. — Nomialized thermal continuum as a function of r/R in the case of the sinusoidal temper- 
ature profile with l/R = 1 for (a) k^R = 10"^ and (b) LR = 10"^ The vertical dotted lines denote 
the boundaries between the different parametrizations of the radiation function. Each zone is in- 
dicated by a number from 1 to 5, that corresponds to the different regimes of Hildner's piecewise 
fit (Table [U). The horizontal dotted line represents 5 = 0. The different line styles in zone 1 corre- 
spond to the Prominence- 1 . 1 (solid), Prominence- 1.2 (dashed), and Prominence- 1.3 (dot-dashed) 
parametrizations. 

The roots of C, = are then = — r-^k-^^ which correspond to the ideal slow continuum, while 
the thermal continuum disappears and becomes the trivial solution 5 = 0. 



3.1. Approximate growth rates and stability criterion 



An approximation to the thermal continuum growth rate can be obtained by neglecting the 
terms with and in Equation ([261 ). Then, the equation = gives the approximate solution 



S ~ Kuk^ + OJt - OJn 

7 \Po 



(28) 



Equation (IzSl) is equivalent to Equation (4.12) of iVan der Linden & GoossensI (119911). It is w orth 
mentioning that Equation (l28l) is similar to the expression obtained by ICarbonell et al.l (|2009|) and 
Solerl (|2010|) for the imaginary part of the frequency of propagating thermal waves in a flowing 
medium. 

The sign of the constant term in Equation (|26T ) provides us with the instability criterion of 
the thermal continuum, i.e., the combination of parameters that causes s > 0. When we use the 
definitions of cOp and coj (Equations (l22l) and (l23l) ) and take into account that the cooling function 
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at equilibrium is L(po, Tq) = when k± = 0, the instability criterion is 



2 (dL podL\ 



(29) 



which turns out to be the same instability criterion given by Field (|1965|) in his Equation (25a), 
although Field's criterion was derived for thermal modes in a homogeneous medium and our equi- 
librium quantities depend on r. Therefore, the stability of the continuum also depends on r, and 
Equation (1291 ) must be computed for all r to assess the absolute instability of the c ontinuum. In 
additio n, magnetic diffusion may play a role for stability. This issue was addressed by llreland et al. 
(|1992l) . who showed that, when magnetic diffusion is included, the contribution of Ohm's heating 
to the equilibrium energy balance has to be taken into account. However, since the magnetic field 
is homogeneous in our case. Ohm's heating plays no role in the energy balance of our equilibrium 
(Equation This fact points out that, in our present application. Equation (|29l ) remains valid 
even when magnetic diffusion is present, and the thermal continuum growth rates are not modified 
by diffusion. 

Equation (1291 ) can be rewritten in terms of parameters x* and a of the cooling function, namely 



(30) 



In the absence of parallel thermal conduction, i.e., k\\ = 0, the instability criterion is satisfied 
for a < 1 . According to the different radiative regimes of Table [1] the condition or < 1 takes 
place for Tq > 80,000 K in Hildner's parametrization, meaning that the thermal continuum is 
unstable in the PCTR-3, PCTR-4, and Corona-5 regimes. Considering the Klimchuk-Raymond 
fit, instability is present for Tp > 93,325 K. This result agrees with the instabilities found by 
Van der Linden & Goossensl(|l99lh in the range of coronal temperatures studied by these authors. 
When parallel thermal conduction is present, the unstable part of the continuum is stabilized by a 
critical value of the longitudinal wavenumber, namely k* given by 



K = 



1 



a 



2 * T*! 



1/2 



(31) 



Note that k* is also a function of r and Equation (I3TI) only applies when a < 1. The absolute stabil- 
ity of the continuum is guaranteed for k^ larger than the maximum value of k* in the equilibrium. 



3.2. Properties of the continuum 



Figure |2] shows the thermal continuum computed from Equation (|28]) for our equilibrium con- 
figuration in the case of the sinusoidal temperature profile with l/R = 1, and considering Hildner's 
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Fig. 3. — Maximum growth rate of the normalized thermal continuum as a function of kji. The dif- 
ferent line styles represent the result for the linear (solid), sinusoidal (dashed), and Gaussian (dot- 
dashed) temperature profiles with l/R = 1. The vertical dotted line corresponds to the maximum 
value of the dimensionless critical wavenumber k*,R for the linear temperature profile. Hildner's 
parametrization for the radiative loss function has been used. 
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Fig. 4. — Same as Figure |2 a) but using Klimchuk-Raymond's fit for the radiative loss function. 

parametrization for the radiative loss function. The growth rate has been plotted in dimensionless 
form computed as 5 = sR/c^p, where s is the dimensionless growth rate and Cgp is the sound speed 
of the central, homogeneous part of the thread. At first sight, we notice the jumps of the growth rate 
at the boundaries of the zones where different parametrizations of the cooling function are used. 
According to Hildner's fit (see Tabled]), the cooling function L(p, T) is continuous at the bound- 
aries between the zones with different radiation regimes, but the derivatives of L (p, T) with respect 
to density and temperature are discontinuous. As Equation (|28]) depends on these derivatives, the 
thermal continuum has the apparent form of five separated continuous spectra. If the actual cool- 
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Fig. 5. — Temperature perturbation (in arbitrary units) of the four most unstable modes (from 
the left to the right) in the absence of magnetic diffusion. Computations performed considering the 
sinusoidal temperature profile with l/R = 1, k^R = 10 ^ and m = 0. The corresponding normalized 
growth rate, s, is displayed on top of each panel. 



ing function of prominences is known and is used here instead of the present parametrization, the 
jumps would be replaced by continuous but very abrupt variations of the growth rate at the bound- 
aries between the different regimes. We must point out tha t the jumps in the thermal s pectrum do 
not affect the stability of the solutions (see Figs. 1 and 3 of IVan der Linden et al. Il99l l. where this 
issue is also comm ented). An alternative radiative cooling function computed numerically (e.g., 



Schure et al.ll2010|) instead of the present piecewise parametrization might be considered in future 



applications. 

In Figure [2t a), corresponding to k^R = 10~^, we see that 5 > in zone 3 and part of zone 4, 
meaning that the continuum is unstable in these regions. This is consistent with the result from 
the instability criterion (Equation (I3U])). For the particular longitudinal wavenumber considered 
in Figure |2ta) parallel thermal conduction can suppress the instability in zone 5 and part of zone 
4. If the longitudinal wavenumber is increased to k^R = 10"' (see Fig. |2tb)), the instability is 
also suppressed in zone 4 and part of zone 3. To completely suppress the instability in the whole 
continuum, we must increase the longitudinal wavenumber until the maximum value of k*R, given 
by Equation (|3TI) is reached. Figure |3] displays the maximum growth rate of the spectrum as a 
function of hR. As expected, the maximum growth rate decreases as k,R increases, until absolute 
stability of the thermal continuum is achieved for a critical longitudinal wavenumber. This critical 
wavenumber is in perfect agreement with Equation (1311 . The result for the three different tempera- 
ture profiles does not show significant differences. We must mention that the various parametriza- 
tions of the radiative regime Prominence- 1 , which aim to represent different optical thicknesses 
of the cool prominence plasma, are not relevant for the instability of the continuum. The thermal 
spectrum in the cool part of the equilibrium, namely zone 1, is always stable independently of the 
considered parametrization (indicated by different line styles in Fig. |2]) 
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In Figure |4] we have plotted the thermal continuum growth rate for hR = 10"^ using the 
Klimchuk-Raymond radiative loss function. By comparing Figures |2t a) and ID we see that slightly 
larger values of the growth rate are obtained using Klimchuk-Raymond's fit with respect to the 
values for Hildner's fit. However, the qualitative behavior of the continuum is similar is both 
cases. 



4, PROPERTIES OF THE UNSTABLE THERMAL MODES 

4.1. Effect of perpendicular thermal conduction 

Here, we study how the unstable part of the thermal continuum is modified when perpendic- 
ular thermal conduction and magnetic diffusion are taken into account. Unless otherwise stated, 
Hildner's parametrization for the radiative loss function is used in all the following computations. 
First, we set 77 = and focus our investigation on the effect of perpendicular thermal conduction. 
We would like to stress that the real physical value of k i given by Equ ation ^ is used in the 



following computations. As stated by I Van der Linden & GoossensI (Il99l|) . the thermal continuum 



is replaced by a set of discrete modes when K_i_ 0. The eigenfunctions of these solutions dis- 
play large variations in a region surrounding the position of the thermal continuum singularity for 
a:^ = 0. Figure |5] displays the temperature perturbation, T^, of the four most unstable modes of 
our equilibrium in the case of the sinusoidal temperature profile with l/R = 1, hR = 10"' and 
m = 0. For simplicity. Figure |5] shows the temperature perturbation only, because Ti is between 
one and two orders of magnitude larger than the other perturbations. This means that the tempera- 
ture perturbation is the dominant disturbance related to the thermal modes, although these solutions 
produce also velocity and magnetic field perturbations. Figure |5] focuses on the r egion where the 



eigenf unctions show significant variations, i.e., the conductive layer described by IVan der Linden 



d 19931) . whereas their amplitude outside the range plotted in Figure |5]is negligible. When Figure |5] 



is repeated for the linear and Gaussian temperature profiles, we find that the eigenfunctions are 
shifted toward smaller r/R for the linear profile, and to larger r/R for the Gaussian profile. How- 
ever, the form of the perturbations is very similar to those displayed in Figure|5]and, for simplicity, 
we do not plot the perturbations again. 

The order of the solutions can be easily identified by the number of extrema (maxima and 
minima) of their temperature perturbation. Thus, the most unstable mode has one maximum only, 
the second most unstable mode has one maximum and one minimum, and so on. The position 
of the largest extremum is shifted toward larger values of r/R as the order of the mode increases. 
In addition, the growth rate (indicated at the top of the panels of Fig. [S]) decreases with the order 
of the mode. These two results are represented together in Figure |6l which displays the growth 
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rate of the 20 most unstable modes as a function of the position of their largest extremum. In 
comparison with the thermal continuum, slightly smaller values of the growth rate are obtained 
for the quasi-continuum modes. The displacement of the largest extremum of the discrete modes 
toward larger values of r/R as their growth rate decreases is consistent with the behavior of the 
thermal continuum. From Figured we also see that the solutions are closer to each other as their 
order increases. This result can be explained by taking into account that the perpendicular thermal 
conductivity (Equation ([7])) is ~ plT~^'^. As we move toward larger r/R, Tq increases and po 
decreases in the equilibrium, meaning that k^^ gets smaller. Thus, the characteristic perpendicular 
spatial-scale decreases, causing the modes to cluster as r/R increases. 

To shed more light on the behavior of the eigenfunctions with k^, we perform the substitution 
Afx Akj^ in the basic equations, with A an enhancing factor. We can artificially increase the value 
of /Cx by means of A and assess its influence on the eigenfunctions. We restrict ourselves to the most 
unstable mode. Figure l^a) shows the temperature eigenfunction of the most unstable mode for 
different values of A. As A grows, the temperature eigenfunction becomes broader. To quantify this 
effect, we compute the width 6 of the maximum of the temperature perturbation measured at its half 
hei ght. The parameter 6 is related to the thickness of the conductive layer. Following the m ethod 
by ISakurai et al.l (|1991|) originally used for Alfven and slow resonances, IVan der Linden il993.) 

1 /3 

obtained that the th ickness of the conducti ve layer is proportional to a:/ . In order to compare the 
analytical result of Van der LindenI (|l993|) with our numerical computations. Figure Ul^) shows 
6/R versus the enhancing parameter A. We see that the dependence of 6/R with A is consistent 
with a scaling law of A^^^, pointi ng out that our numeric al code works properly and recovers the 
behavior analytically predicted by IVan der LindenI (|l993|) . 



We have also computed the eigenfunctions for other values of the azimuthal wavenumber 
m. We do not have obtained variations of the spectrum of quasi-continuum modes for different 
values of m. It is worth noting that for m = 0, the perturbations and are decoupled from 
the remaining perturbations, while for m all perturbations are coupled. The temperature 
perturbation is independent of m. The effect of m on the growth rate is studied in Section |5] 



4.2. Effect of magnetic diffusion 



Here, we firs t consider magnetic diffus ion and neglect perpendicular thermal conduction, i.e., 
A'x = and rj i^O. Ilreland et al.l(ll992lll998|) studied the effect of magnetic diffusion on the thermal 
continuum. For their particular equilibrium, these authors studied two limit cases of the relative 
values of the thermal and Alfven time scales. W hen the thermal time scale is much longer than 
the Alfven time scale. Ilreland et al.l (|1992lll998[) found that the effect of diffusion is to replace the 
thermal continuum by a dense set of discrete modes, whose growth rate is displaced with respect 



-19- 



0.30 




0.15 - 



0.1 OL . . . . : . . . . : . . . . J 

0.675 0.680 0.685 0.690 

r/R 



Fig. 6. — Normalized growth rate of the 20 most unstable modes (symbol x) as a function of the 
position of their largest maximum or minimum. The solid line represents the thermal continuum 
in the case = 0. Results are computed for rj = 0, kji = 10 \ and m = 0. The sinusoidal 
temperature profile with l/R = 1 has been used. 
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Fig. 7. — (a) Temperature perturbation (in arbitrary units) of the most unstable mode. The different 
line styles correspond to different values of the perpendicular conductivity enhancing factor: X = \ 
(solid), /I = 10 (dotted), and /I = 100 (dashed), (b) Width of the maximum of the normalized 
temperature perturbation, 5IR, corresponding to the most unstable mode (symbol x) versus the 
perpendicular conductivity enhancing factor X. The dotted line corresponds to the dependence 
6IR ~ A^^^. Results are computed for i] = 0, k^R = 10~\ and m = 0. The sinusoidal temperature 
profile with l/R= I has been used. 
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Fig. 8. — (a) Ratio of the Alfven time scale, ta, to the thermal continuum time scale, Tj- The 
vertical dotted lines denote the boundaries between the different parametrizations of the radiation 
function according to Table [T] (b) Normalized growth rate of the most unstable mode versus the 
diffusivity enhancing factor e. The sinusoidal temperature profile with IjR = 1 has been used. 



to the values of the thermal continuum. The displacement of the growt h rate is independent of th e 
value of the diffusivity. This is the case of the cool profile studied by Ireland et al. ( 1992 . 1998|) . 
On the other hand, if the thermal time scale is much shorter than the Alfven time scale, the thermal 
continuum is not affected by di ffusion and the discrete m odes are absent. This other possibility 
corresponds to the hot profile of Ireland et al. ( 1992 , 1998() . 



To assess the effect of magnetic diffusion on the thermal continuum of our equilibrium, let 
us compare the Alfven time scale, ta = Ljvf^, with the thermal time scale, tj = l/\s\, where s 
is here the thermal continuum growth rate and L is a typical length scale. We relate L with the 
longitudinal wavelength of the perturbation, namely L = In/h. This is done in Figure[8]^a), which 
displays the ratio ta/tt versus r/R with hR = 10"'. We obtain that in the most unstable part of 
the continuum, i.e., the beginning of zone 3, Tr and ta are o f the same or der. Therefore, we are 
in a situation between the limit cases studied by Ireland et al.l (119921 119981) . We have numerically 
checked that the growth rates are not modified by the presence of magnetic diffusion. In addition, 
by setting rj and k± = in the equations implemented in the PDE2D code, we have found no 
discrete t hermal solutions. Therefo re, the behavior of our equilibrium is similar to that of the hot 
profile of llreland et al.l (|l992lll998b although in our case Tt and ta are of the same order. 



A different issue is to determine how the previously described modes for /c^ ^ and rj = 
(Sec 14.11) are affected when magnetic diffusion is included. Hence, we consider the general 
case /Cx ^ and rj 0. As was done for the perpendicular thermal conduction, we perform the 
substitution rj ^ erjin the basic equations, with e an enhancing factor. We can artificially increase 
the value of 77 by means of e and assess its influence. As before, we restrict ourselves to the most 
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unstable mode, whose growth rate versus e is displayed in Figure [8tb). For e < 50 the growth rate 
is constant and coincides with the value for 77 = 0, while the growth rate increases for e > 50. For 
such large values of e, magnetic diffusion govems the behavior of the solutions. When the actual 
value of J] is considered, i.e., 6=1, the growth rate is approximately the same as that obtained 
in the absence of magnetic diffusion. Thus, we need much larger, unrealistic values of rj for the 
growth rate to be affected by diffusion. On the other hand. Figure |9] shows the evolution of the 
temperature perturbation of the most unstable mode as e increases. Again, the result for e = 1 is 
almost indistinguishable from the eigenfunction for 77 = (compare Figs. |5la) andjOfa)). When e 
is increased and so very large, unrealistic values of rj are considered, the eigenfunction is affected 
by an additional modulation caused by magnetic diffusion and develops smaller spatial scales as e 
is increased. On the basis of these results, we conclude that, for realistic values of the diffusivity, 
magnetic diffusion is irrelevant for both the growth rate and the eigenfunctions of the most unstable 
modes. 
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Fig. 9. — Temperature perturbation (in arbitrary units) of the most unstable mode for a diffusivity 
enhancing factor of (a) 6 = 1, (b) e = 50, (c) e = 70, and (d) e = 100. Computations per- 
formed considering the sinusoidal temperature profile with l/R = 1, hR = 10 \ and m = 0. The 
corresponding normalized growth rate, s, is displayed on top of each panel. 



4.3. Comparison between the Hildner and Klimchuk- Raymond radiative loss functions 

We have performed some test computations using Klimchuk-Raymond's fit for the radiative 
loss function in order to compare with the results obtained using Hildner's fit (discussed in the 
previous Subsections). Regarding the form of the temperature perturbations, we find no signifi- 
cant differences. The maximum of the temperature perturbation is shifted toward slightly larger 
values of r for Klimchuk-Raymond's fit. For the fundamental mode, the maximum of takes 
place at r/R x 0.675 using Hildner's fit, and at r/R x 0.69 using Klimchuk-Raymond's fit for 
kyR = 10"\ l/R = 1, and m = 0. On the other hand, the growth rate is larger for Klimchuk- 
Raymond's parametrization. This is in agreement with the behavior of the thermal continuum (see 
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Fig. 10. — Normalized growth rate of the most unstable mode as a function of (a) the azimuthal 
wavenumber, m, for kji = 10"\ and (b) the dimensionless longitudinal wavenumber, kM, for 
m = 0. The results of both panels correspond to the sinusoidal temperature profile with l/R = 1. 




Fig. 11. — Normalized growth rate of the most unstable mode versus the thickness of the tran- 
sitional layer, l/R. The different line styles in panel, almost superimposed, represent the linear 
(solid), sinusoidal (dashed), and Gaussian (dot-dashed) temperature profiles. All computations are 
performed with k^R = 10^ and m = 0. 

Figs. 121; a) and HI). For the same set of parameters as before, the growth rate of the most unstable 
mode is 5 « 0.1994 using Hildner's fit, and s ~ 0.4370 using Klimchuk-Raymond's fit. Apart 
from this difference in the growth rate, the results for both Hildner's and Klimchuk-Raymond's 
parametrizations are equivalent. 
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5. PARAMETRIC STUDY OF THE GROWTH RATES 

In this Section, we study the growth rate of the most unstable mode as a function of the 
dimensionless longitudinal wavenumber, k^R, the azimuthal wavenumber, m, and the thickness of 
the transitional layer, IjR. Hildner's fit for the radiative loss function is used in all cases. 

First, we assess the efl'ect of the azimuthal wavenumber, m, (see Fig. [TUta)). Note that only 
positive integer values of m are considered in Figure [ISa) as the results are equivalent when 
negative values are used. As for the thermal continuum, the value of m is irrelevant for the growth 
rate of the discrete thermal modes. The almost negligible variation of the growth rate with m shown 
in Figure [Tola) might not be physical and might be attributed to the accuracy of the numerical code. 
On the other hand, the dependence on the longitudinal wavenumber is displayed in Figure [TOl b). 
The maximum of the growth rate takes place for hR « 0.4, while s decreases as kji is increased 
or reduced. For kJi > 0.4, the growth rate decreases and the thermal mode is stabilized by the 
eff'ect of thermal conduction parallel to magnetic field lines, as commented in Section [3l The 
critical wavenumber for the stabilization is consistent with the expression given in Equation (I3TI ) 
for the stabilization of the thermal continuum. The growth rate also decreases with respect to 
the maximum value when k,R < 0.4. This behavior is different from the behavior of the thermal 
continuum (compare Figs . |3] and [TOf b)) . In addition, we find that for k,R ^ a rich collection of 
couplings between the diff"erent thermal modes takes place. This complex system of couplings is 
not visible on the scale of Figure (TOfb). We do not explore these couplings in detail here because 
they take place for very small kJi and far from the maximum value of the growth rate. 

Finally, the dependence on the thickness of the transitional layer, l/R, is shown on Figure [TTI 
In a real prominence thread, the thickness of the transitional layer is probably determined by the 
coefficient of thermal conduction perpendicular to the magnetic field. In fully ionized coronal 
plasmas, the perpendicular thermal conductivity is very small. However, for cool prominence tem- 
peratures the plasma is only partially ionized and the perpendicular thermal conductivity in creases 



by several orders of magnitude due to the effect of th ermal conduction by neutrals (see, e.g.. Parker 



19531 : llbanez & Mendozalll990l : iForteza et al.ll2008[) . This effect might be important in the inner 



part of the PCTR. Unf ortunately, since the ionization degree in prominence threads is unknown 



(see recent estimations looutteb roze & Labrosse 2009 ). the thickness of the transitional layer can 



not be determined with confidence. For this reason, we take l/R as a free parameter. The growth 
rate is almost constant for l/R > 0.5, whereas s decreases for l/R < 0.5. When l/R — > 0, S — > 
because the transitional layer is absent and the thermal continuum disappears. As in previous com- 
putations, the results for the different temperature profiles are almost identical, hence the form 
of the temperature profile is not relevant for the growth rate. This is an important result because 
it means that our present conclusions apply for all possible temperature profiles in the particular 
PCTR of prominence threads. 
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6. IMPLICATIONS FOR THREAD LIFETIMES 

In this paper we have shown that prominence threads, modeled as cool and dense magnetic 
flux tubes embedded in a much hotter and less dense coronal environment, are not stable because 
unstable thermal mo des are presen t in their particular PCTR. This statement qualitatively agrees 
with the results from Karpen et alj ([l989), who studied thermal instabilities in a two-dimensional 



spatially random magnetic field distribution, and concluded that cool condensati ons cannot be in 



static equilibrium with a hot exterior. However, in the work by lKarpen et alj (|1989|) the surrounding 
of the condensations simply continue to cool down because these authors did not consider the 
eff"ect of thermal conduction perpendicular to the magnetic field. The unstable discrete modes 
present in our model appear due to the eff'ect of thermal conduction perpendicular to the magnetic 
field. Our results for the growth rate of the unstable thermal modes have direct implications for the 
stability and lifetime of prominence threads. In our analysis, we have restricted ourselves to the 
linear stage of the thermal instability. This stage represents the initial phase in which the unstable 
thermal mode grows after its excitation. The combined influence of the diff"erent unstable modes 
will afl'ect the subsequent nonlinear evolution of the plasma within the PCTR, because many of 
the unstable modes have very similar growth rates and the evolution will not be dominated by a 
single mode. Significant changes in the equilibrium of the threads would take place during the 
subsequent nonlinear evolution of the thermal instability, because of the large temperature and 
density gradients generated, which may be detected by the observations. The growth rate of the 
linear phase provides us with a reasonable estimation of the typical time scale on which the efl'ect 
of the thermal instability in the PCTR of the threads would be observable and may afi'ect the thread 
dynamics. 

In the analysis of the previous Sections, we have seen that the most relevant parameter afl^ect- 
ing the instability growth rate is the longitudinal wavenumber, since it determines the efficiency of 
parallel thermal conduction as stabilizing mechanism. If unstable thermal modes are excited in a 
prominence thread by an arbitrary disturbance, the instability would be dominated by the most un- 
stable solution. We find that in our model the maximum of the growth takes place for k^R x 0.4, and 
corresponds to a dimensionless growth rate of 5 0.255 using Hildner's radiative loss function, 
and s ~ 0.467 using Klimchuk-Raymond's radiative loss function. The corresponding instabil- 
ity time scale is = l/s, with s the dimensional growth rate computed from its dimensionless 
value as 5 = sCsp/R. In Figure [I2l we plot tt as a function of the temperature at the center of the 
thread, Tp, for different values of the thread radius, R. We find that, for realistic values of Tp and 
R, the instability time scale is of the order of a few minutes. We obtain that tt decreases as Tp 
gets higher and increases as R gets larger, meaning that cool and wide threads are more thermally 
stable than hot and thin threads. We also see that Klimchuk-Raymond's parametrization produces 
smaller instability time scales than with Hildner's parametrization. The values of tj indicated in 
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Figure[T2]are of the same order of magnitude as the t ypically observed lasting ti r ne of the threads 
in sequences of Ha images of solar filaments (e.g., Linl 120041: Lin et al.l l2005l |2009|) . This re- 
sult implies that the thermal instability may play a relevant role for the dynamics and stability of 
prominence and filament fine structures, because the time scale related to the thermal instability is 
consistent with the observed time scale on which threads seem to suffer important changes in their 
equilibrium. 
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Fig. 12. — Thermal instability time scale, tj-, computed from the maximum growth of the most 
unstable mode versus the temperature of the thread core, Tp, using (a) Hildner's and (b) Klimchuk- 
Raymond's parametrization for the radiative loss function. The different line styles in both panels 
represent different values of the thread radius: R = 100 km (solid), 200 km (dotted), 400 km 
(dashed), and 800 km (dot-dashed). 



7. DISCUSSION AND SUMMARY 



In this paper we have studied the properties of quasi-continuum thermal modes in a promi- 
nence thread model which is tr ansversely inhomogeneous. We h ave followed the method of 
Van der Linden et al.l (|l99l|) and IVan der Linden & GoossensI (|l99l|) to investigate, first, the sta- 
bility of the thermal continuum and, later, the effec t of cross-field therma l condu ctio n and magnetic 

diffusi on. We have recovered the general results by IVan der Linden et al.l (|l99lh and lVan der Linden & Goossens 
(|199l h in all cases studied in the present paper. In particular, we have found that the thermal con- 
tinuum in prominence thread s is un stable for PCTR temperatures. In agreement with the results of 
Van der Linden & GoossensI (|l99l|) . the effect of cross-field thermal conduction is to replace the 
thermal continuum by discrete modes which retain the unstable characte r of the continuum. On th e 
contrary, the role of magnetic diffusion is negligible in our model (see llreland et al.lll992l Il998|) . 
We have obtained the important result that the instability growth rate of the most unstable mode is 
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independent of the fomi of the temperature profile within the PCTR of the thread, and the insta- 
bility time scal e is consistent with the observed lifetime of the threads in Ha observations of solar 
fil aments (e.g.jLin et al.ll2008L 120091) . Co nsidering our present r esults along with those obtained 
by lVan der Linden & Goossensl(|l991|) and lVan der Lindenl(|l993|) . we conclude that unstable ther- 
mal modes may play a relevant role for both the formation of the prominence fine structure and the 
subsequent instability of the thin prominence threads. 

Here we have considered a fully ionized prominence. However, due to the low temperature 
in the dense core of the thread s, the cool prominence plasma is expected to be partially ionized 
( Gouttebroze & Labrosse 200 91). If partial ionization is taken into account, several additional ef- 
fects have to be considered, namely thermal conduction by neutrals, ambipolar (or C owling's) dif- 
fusion caused by ion-neu tral collisions, and a modified radiation function (see, e.g., Forteza et "aD 
20081 : ISoler et al.ll2010al) . The role of both effects depends strongly on the plasma ionization de- 
gree. It requires complicated computations of NLTE radiative transfer and statistical equilibrium of 
atomic level populations to ob tain the precise profile of the ionization degree in prom inence threads 
(see extensive details in, e.g. jGouttebroze & Labrossell2009l : iLabrosse et al.ll2010l) . Nevertheless, 
some relevant conclusions can be obtained if one reasonably assumes that the ionization degree in 
prominence plasmas is mainly determined by the temperature and, to a lesser extend, by the den- 
sity. While the ionization degree of the core of the thread is unknown, it is realistic to assume that 
the condition of full ionization o ccurs in the PCTR for temperatu res higher than a c ritical temper- 
ature, namely T* . The results by iGouttebroze & Labrossel (|2009|) and lSchure et al.l (120101) suggest 
realistic values for the critical ionizing temperature for the hydrogen plasma in prominences of 
T* X 2x 10^ K and T* ^ 2.5 x lO'^ K, respectively, so only the inner part of the transitional layer 
would be partially ionized. According the instability criterion (Equation (|30] )). the unstable part 
of the thermal continuum takes place at temperatures higher than 8 x 10*^ K, which is higher than 
the critical ionizing temperature. This suggests that the ionization degree of the cool part of the 
fine structure is probably not relevant for the unstable quasi-continuum modes, because they are 
mainly confined where hydrogen may be already fully ionized. However, a detailed investigation 
of this issue is needed for more robust conclusions. 



The present investigation is a first step and should be improved in the future by taking into 
account relevant effects that have not been included in our model. Among these effects, the role of 
density and magnetic longitudinal inhomogeneity and the effect of magnetic twist will be assessed 
in forthcoming works. We have also neglected the effect of mass flows. Pr evious studies on 
linear thermal modes in flowing homogeneous prominence plasmas (see, e.g.. ISoler et al.l 12008 : 
Carbonell et al.ll2009l 1201 Oe iBarcelo et al.ll201 1|) suggests that flows along magnetic field lines do 
not affect the linear growth rate of the thermal instability. 
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